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Abstract 

In this letter we continue the calculation of master integrals for massless three-loop form 
factors by giving analytical results for those integrals which are relevant for the fermionic 
contributions proportional to Np, Np-N, and Np/N. Working in dimensional regularisation, 
we express one of the integrals in a closed form which is exact to all orders in e, containing 
T-functions and hypergeometric functions of unit argument. In all other cases we derive 
multiple Mcllin-Barnes representations from which the coefficients of the Laurent expansion 
in e are extracted in an analytical form. To obtain the finite part of the three-loop quark and 
gluon form factors, all coefficients through transcendentality six in the Riemann ^-function 
have to be included. 



1 Introduction 



The quark form factor 7* — > and gluon form factor H ^ gg (effective coupling) are the 
simplest processes containing infrared divergences at higher orders in massless quantum field 
theory, and therefore appear in a large variety of physical applications. They can for instance 
be used to predict the infrared pole structure of multi-leg amplitudes at a given order [1-3] and 
to extract resummation coefficients [4,5], and they are needed for the purely virtual corrections 
to a number of collider reactions (Drell-Yan process, Higgs production and decay, DIS). 

At the two-loop level, corrections to the massless quark [6] and gluon [7, 8] form factors were 
computed in dimensional regularisation with D = 4 — 2e to order e° and subsequently extended 
to all orders in e in Ref. [9]. Two-loop corrections to this order were also obtained for massive 
quarks [10]. The three-loop form factors to order e _1 (and e° for fermion loop contributions) 
were computed in [5,11]. One of the main motivations for obtaining analytical results for the 
form factors is the search for a deeper underlying structure of the coefficients, as proposed in 
Ref. [12] for planar box amplitudes. 

In order to calculate the quark and gluon form factors at higher orders in perturbation theory, 
the amplitudes are reduced, by means of algebraic reduction procedures [13-16], to a small set 
of master integrals. At the three-loop level, the master integrals for massless form factors were 
identified in Ref. [17] and results for certain subsets are available in the literature [13,17-19]. 
Among the three-loop master integrals, the genuine three-loop vertex functions are the most 
challenging ones from a computational point of view. They correspond to two-particle cuts of 
the master integrals for massless four- loop off-shell propagator integrals [20], which have been 
used in the calculation of the scalar i?-ratio [21]. The derivation of the three-loop vertex integrals 
is of comparable complexity to the four-loop propagator integrals. 

Working in dimensional regularisation and expanding the master integrals in a Laurent series 
in e, the finite part of the three-loop form factors requires the extraction of all coefficients through 
transcendentality six, i.e. coefficients containing terms up to 7r 6 or (3. Note that the power of 
e coming with coefficients of transcendentality six in the Laurent expansion is not always the 
same in the different master integrals: Transcendentality six can appear in the coefficients of 
the e -, e- or e 2 -terms in the Laurent series. If it appears in the e fc -term, this indicates that a 
prefactor ~ l/e fc will come from the reduction to master integrals, such that an expansion up 
to transcendentality six of the master integrals will always be required. 

Those genuine three-loop vertex functions which contain one-loop or two-loop propagator 
insertions were already given in Ref. [17]. The purpose of the present letter is to extend this 
calculation to all three-loop master integrals which have less than nine propagators. Each 
topology contains only one master integral, which is chosen to be the scalar integral, with no 
loop momenta in the numerator and with all propagators raised to unit power. It turns out 
that this subset of three-loop master integrals is sufficient in order to obtain the aforementioned 
fermion loop contributions within a Feynman diagrammatic approach [22]. At this point we 
would like to point out an error in Ref. [17], namely the basis of three-loop master integrals 
given there is too large, since certain two-particle cuts of four-loop propagator topologies [20] 
yield topologically identical three-loop vertex topologies. Consequently = Ag^ = ^8 an d 
Ag^ = Aq^. The corrected set of three- loop vertex integrals is given in Fig. [TJ 

This letter is organised as follows. Computational methods to obtain analytical and nu- 
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merical results of the three-loop vertex integrals with up to eight propagators are described in 
Section [21 and the analytical results for them are listed in Section [3j Section 2] contains our 
conclusions and an outlook. 

2 Master integrals: Classification and computational methods 

Vertex integrals with one off-shell and two on-shell legs and massless propagators depend only 
on one kinematic scale: the mass q 2 of the off-shell leg. The dependence on this scale is given 
by the mass dimension of the integral, such that the coefficients of the Laurent expansion in 
the dimensional regularisation parameter e are real constants which are in general of increasing 
transcendentality in the Riemann (^-function, where the degree of transcendentality (DT) is 
defined by 

DT(r) = for rational r 
DT(ir k ) = DT(((k)) = k 
DT(x ■ y) = DT{x) + DT(y) . (1) 

At the three-loop level the quark form factor depends - like the process e + e~ — > 3 jets at NNLO 
- on the following seven colour structures [22,23] 

N 2 , N°, 1/N 2 , N F -N, Np/N , N 2 , N F>y , (2) 

where the last colour factor stems from topologies in which the external gauge boson couples to 
a closed fermion loop. The three terms containing Np are referred to as fermionic corrections. 
They have been derived in Refs. [5, 11] from the behaviour of the three-loop deep inelastic 
coefficient functions [24] . In the more conventional approach of computing multi-loop Feynman 
amplitudes the form factors are - after an algebraic reduction procedure [13-16] - expressed in 
terms of a small set of master integrals. It turns out [22] that the master integrals in Fig. [1] with 
at most eight propagators are sufficient in order to obtain the fermionic corrections to the form 
factor. The purpose of this letter is therefore to evaluate these master integrals. Those master 
integrals in Fig. Q] that contain single or multiple bubble insertions have already been computed 
in Ref. [17], the remaining ones with up to eight propagators - i.e. diagrams Aq^, ^7,3, ^4.7,4, 
A75, and A$ - are subject of the present work. Working in dimensional regularisation with 
D = 4 — 2e, we give one of the diagrams (^7,4) in a closed form which is exact to all orders 
in e, containing T-functions and hypergeometric functions of unit argument. In all other cases 
we derive multiple - twofold to fourfold - Mellin-Barnes representations [25-28] from which the 
coefficients of the Laurent expansion in e are obtained in an analytic form. As explained above, 
all coefficients through transcendentality six in the Riemann ^-function have to be included to 
obtain the finite part of the three-loop form factor. 

For many practical applications, and to verify the analytical results, it is sufficient to know 
the numerical values of the coefficients in the Laurent expansion of the master integrals to some 
finite order. There are several techniques to obtain numerical values for the coefficients, one 
of them being the sector decomposition method, which is described in detail in Refs. [29,30]. 
Using this technique, the Laurent expansions of all master integrals relevant to the three-loop 
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form factors can be computed to, in principle, any desired order. In practice there are of course 
limitations, from CPU time for the numerical evaluation and from memory for the algebraic part 
of the sector decomposition procedure. The eight propagator graph A$ is the most complex one 
from the sector decomposition point of view, not only due to the high number of propagators, 
but also because it exhibits spurious linear divergences at intermediate stages, which render the 
subtractions and thus the functions to be integrated more complicated. The computing time for 
A$ up to order e for a numerical precision of 0.1% is of the order of 4 hours on a 3.0 GHz PC. 
For a precision of 1% the evaluation is more than 10 times faster. 

Another method of doing numerical cross checks proceeds along the lines of deriving Mellin- 
Barnes representations of loop integrals by means of the package AMBRE [31] and subsequently 
performing the analytical continuation and numerical evaluation of the obtained expressions 
with the package MB [32]. This procedure allowed us to check most of the coefficients at the 
sub-permille level. 



3 Results 



In this section we list the results we obtained for the three-loop master integrals necessary for 
the fermionic corrections to the three-loop quark form factor, which are the diagrams Aq^, Aj^, 
^7,45 ^7,5 an d As i n Figure [H All other diagrams with up to eight propagators possess so-called 
bubble insertions and have already been given in Ref. [17]. 



Diagram A 



6,2 



The first diagram to be considered is ^6,2- hi Ref. [33] a representation of this diagram in 
terms of a one-dimensional integral over hyper geometric functions was given. Here we pursue 
a different strategy and derive a twofold Mellin-Barnes representation [26-28] from which the 
coefficients of the Laurent series expansion about e = can be computed. We start with 



A 



6,2 



d D k r d D i r d D 



(2ir) D J (2ir) D J (2ir) D (k + Pl f {k + I - p 2 f I 2 r 2 (r - k) 2 (r - k - I) 



2 " 



(3) 



and assume here and in the following that all propagators contain an infinitesimal +ir]. We then 
derive the following expression that contains a triple integral over a Meijer-G function [34,35], 



A, 



6.2 



i rj 



-3e r 3 (l-e)r(3e) 
r(l-2e)r(2-4e) 



J dx dy dz x~ e (1 



)- 3e y~ e (1 - y)-' Se z' 2t (1 - z) 



3 c 



-2c 



where 



xG§(xz + y(l- :) 



(Pi +P2? 



{-l + 4e,-l + 4e} , {3e} 
{-l + 3e,-l + 2e,0} , {} 
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(4vr) D/2 T(l - e) 
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We now make use of the contour integral representation of the Meijer-G function [34,35], and 
subsequently decompose the argument by means of a second Mellin-Barnes representation. The 
integrals over x, y, and z can then be done in terms of T-functions. This leads to the following 
twofold Mellin-Barnes representation for A§ 2 , 



A 



6,2 



1 rj 



-3e r 3 (l -e)T(3e) T 2 (l -3e) 



2e)T(2 - 

r(-l + 3e-w 1 ) r(-l + 2e - i0i)r(2 



4e) 



2vri~ 



du> 2 
2vri" 



4 e + r(— iy 2 ) r(w 2 - u>i) 



r(3e- wi)T(2-4e + u;2)r(2-4e + u;i - iy 2 ) 



X T(l - e + w 2 ) r(l - e + wi - ib 2 ) T(l - 2 e + w 2 ) r(l - 2 e + wi 



W2 



(6) 



1 3 

< e < — 

15 20 



(7) 



In the above equation ([6]) the contour integrals in the complex plane are along curves which 
separate left poles of T-functions from right ones, where "left poles" are poles stemming from a 
r(. . . +w) dependence, while "right poles" stem from a T(. . . — w) dependence [28]. The most 
convenient choice for these contours are straight lines parallel to the imaginary axis, i.e. the real 
parts along the curves are constant. According to Refs. [26,27], these real parts, together with 
the parameter e, must be chosen in such a way as to have positive arguments in all occurring 
r-functions in order to separate left and right poles in the desired way. One verifies easily that 

6 1 

Cl = "5 ' 02 = "2 ' 
is an appropriate choice in Eq. ([6]). From the fact that the origin lies within the allowed region 
for e, we conclude that the Mellin-Barnes integration does not produce any poles in e in addition 
to the UV pole that is already present in the prefactor. Therefore the expansion in e commutes 
with the contour integrations. Proceeding in this way, the Mellin-Barnes integrations can be 
done order by order in e. During this procedure, the contours can be closed at infinity to either 
side of the complex plane, and the corresponding residues are then summed with the appropriate 
global sign. Furthermore, the nested sums algorithm [36,37] and the formulas in the Appendix 
of Ref. [28] - Barnes Lemmata and corollaries thereof - prove extremely useful. The final result 
for A e 2 is 



A 



6.2 



1 r\ 



-3 c 



2C 3 7vr 4 / 7tt 4 2tt 2 h \ 

— - I8C3 + -122 C3 + C3 - IOC5 e 

e s 180 \ s 20 3 / 



+ -738 Cs 



427vr 4 



+ 6^ 2 C 3 -90C 5 + 



1637T 6 



+ 76Cf e 2 + 0(e 3 ) 



(8) 



180 ' ' 7560 

In Ref. [33], two more orders of the e-expansion can be found. 

Diagram Aj^ 

We now turn our attention to the integral This integral will be represented, similarly 



5 



to the integral Aq 2 , as a multiple Mellin-Barnes integral: 
d D k f d D l fd D r 



A 



7,3 



(2tt) D J {2tt) u J {2tt) u k' 2 (k + qf (I - k - p 2 f (I - p 2 f (r + if r 2 (r - Pl f 



i r\ 



-1-3 e r 4 (l — e)r(— e) f dw\ f dw2 f dw^ T(—wi) 



r(l-2e)r(l-3e) J 2m J 2iri J 2m T{1 - wx) 

ci -i 00 c 2~ * 00 eg — ioo 

r(-3 e - w 3 ) T(l + 2 e + Wl + w 2 ) T(l + ^ + u; 2 ) T(-2 e - w 2 ) T(-e - Wl ) 



r(l - 3 e - w 3 ) T(2 - 2 e + + w 2 ) 

x r(-u> 3 ) r(e - iui — W2 + ^3) r(i — e + ^2) r(i + w 3 ) r(— e + - tf 3 ) 



(9) 



The contour integrals are again along straight lines in the complex plane parallel to the imaginary 
axis, and as before we must choose the real parts of the integration variables such as to have 
positive arguments in all occurring T-functions. This is achieved by choosing 



Cl = "20 ' 



C-2 



C3 



1 3 

'8 <£< 20- 



(10) 



As it was the case for ^6,2) we have the origin within the allowed region for e and therefore the 
Mellin-Barnes integration does not give rise to any additional poles in e, the only pole of the 
integral being the infrared pole that is already present in the prefactor in Eq. ([9]). We can thus 
again perform the contour integrations order by order in e. Since the leading coefficient turns 
out to have already transcendentality five, we only need to compute the first two terms in the 
expansion. They are given by 



A 



7.3 



l Tj 



-l-3e 



vr 2 C3 
6 



10 6 - 



1 119 7T b 



2160 



(11) 



During the evaluation of this integral we could not proceed until the end by merely applying 
Barnes Lemmata and corollaries thereof, but had to apply auxiliary integral representations of 
hypergeometric functions at intermediate steps. The benefit of this procedure is that it enables 
us to perform all Mellin-Barnes integrations, at the cost of introducing additional parameters 
over which we subsequently have to integrate. However, the integrations over these auxiliary 
parameters can be done in terms of logarithms and (harmonic) poly logarithms. Furthermore, 
we made extensive use of the package HPL [38,39] and of an algorithm based on the nested sums 
approach [36,37]. See Appendix [A] for more details on this point. 

Diagram Aj 4 

The next diagram we consider is At first glance it looks quite difficult since it lacks 

both a bubble insertion and a planar topology. However, it turns out to be simpler than the 
planar diagram ^7,3, and it can even be displayed in a closed form. The main reason for this is 
the fact that at the outer vertices of both outgoing lines only three lines meet, and hence the 
introduction of Feynman parameters allows for the combination of propagators that differ only 
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by a light-like momentum. This property is absent in both A e 2 and A 73 , both of which did not 
reveal a closed form but only a multiple Mellin-Barnes representation. For A-?^, we find 



rdPk r d D i rd D 7 

7,4 J (2tt) d J (2tt) d J~(2^) 



(2ir) u J (2vr) w J (2ir) v k (k - qf (r + 1 - k) I 2 (I — p\) r 2 (r - p 2 Y 

1-3 e 



X 



r(l-e)r(3e) 



2-r 4 (l-e) T 2 (-e) 



4 -F 3 (l,l-e,l-3e,2-6e; 2-3e,2-3e,2-4e; 1) 



(1 -3e) z T(2-4e) 

T(l - 3 e) T(2 - 3 e) T(3 e) T(l + 2 e) 
r(2-e)r(2-6e) 

x 4 F 3 (l,l,l + 2e,2-3e; 2,2,2-e; 1) 

r 2 (i-3e)r(i + 2g)r(i + 3g) 

+ r(2-e)T(2-6e) 

x 4 F 3 (l,l,l + 2e,l + 3e; 2,2,2-e; 1) 



i r\ 



-l-3e 



X 



— + v ^ + 36C3 J7 + l v — 
/ 22^ 4 19 vr 6 



+ 216C 3 -2^ 2 C3 + 46C 5 



+ 



+ 1296 Cs - 12 vr 2 Cs - 282 Cf + 276 Cs U + 0(e 2 ) 



(12) 
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where the expansion in e was done by means of the Mathematica [40] package HypExp [33,41]. 
Diagram A-j^ 

We now consider diagram A-j^ for which we derive the following fourfold Mellin-Barnes 
representation. 

• d D k f d D l f d D r 1 



fd u k f d D l fd u r 
7,5 JJ^JJ^JJ^ 



(2vr) D J (2vr) D J (2vr) D k 2 (k + qf {k + rf (I - p 2 f (r-lf r 2 (k + I + pif 

r 3 (i - e) 



-l-3e 



c-j^+i oo c 2"l"* 00 co+ioo c^ + i oo 



r(l - 2e)T(l -4e) 7 2vni J 2vri J 2vri J 2m 

ci—ioo C2~ioo C3 — ioo c^—ioo 



T(w 4 - wi) r(l + io 3 ) T(-3e - u> 3 ) T(l - 2e + u>i + w 2 - 104) T(-w 3 ) r(-u> 4 ) 

r(i - wi + w s + 104) 



T(l + 3e + w 3 + Wj) T(l + e + wi + w 2 ) T(l + Wi + w 2 ) r(-e - wi) Fj-e - w 2 ) 
T(2 + e + w 1 + w 2 )T(2-2e + w 1 +w 2 ) 

T(l - e + w 2 ) T(l + w 3 ) T(l - e + ti^) r(e - wi - w 2 + w 3 + w 4 ) T(w4 - w 2 ) 



r(i -w 2 + w 3 + w 4 ) 



(13) 



Like in the previous cases, the Mellin-Barnes integral does not generate poles in e, so we can 
therefore interchange the expansion in e with the contour integrations. We choose 



ci 



1 

5 



C2 



1 



C3 



1 
7 



C4 = - n . 



(14) 



As before in the case of we can not proceed until the end by merely applying Barnes 
Lemmata and corollaries thereof, but again have to apply auxiliary integral and series represen- 
tations at intermediate steps, this time even for a larger class of functions than before. Besides 
hypergeometric functions, these are mainly logarithms and (harmonic) polylogarithms as well 
as V'-functions with 



d_ 

dz 

d_ 

dz 



ln[T(z)] , 

^-^(z) for k=l,2,... 



(15) 



The sums and integrals over the auxiliary parameters are then performed by means of the same 
techniques as before. Some details can again be found in Appendix A. The final result for A?fi 
reads 



A 7 , 5 = i S% 



Diagram A$ 



1 77 



-l-3e 



27T 2 Cs + 10 C5 + ( 12^ 2 C3 + 60 C 5 + ^ + 18 C 3 2 ) e + 0(e 2 ) 



(16) 



The last diagram we consider is A$ which can also be displayed as a fourfold Mellin-Barnes 
integral. 



A* = 



d D k f d u l [d u r 



(2tt) d J (2tt) u J (2tt) u {k + Pl f (k + rf {k + r + qf (I — k) z (I + rf I 2 r 2 (I + Pl 



1 



2 7 2 r 2 



— q —irj 



-2-3e T 3 (l - e) r(-l - 3e) / dwi f dw 2 f dw. 



/dw\ f dw 2 f du>3 f du>4 
2^i J 2ri J ~2~^i J 2ri 



r(l + w 3 ) T(l + W4) T(w4 - w 2 ) T(w 3 - w\) T(-W4) T(—w 3 ) T(2 + w\ + w 2 ) 

r(2 + w 3 + W4) r(i + W4- w 2 ) r(i + ^ 3 - wi) 
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r(2 + e + W! + w 2 ) T(l - e + Wl ) T(l - e + w 2 ) T(-l - e - Wl ) - e - w 2 ) 
T(2 - 2e + u-i + w 2 ) T(3 + e + u-i + w 2 ) 



x T(2 + 3e + W3 + u>4) T(l — 2e + wi + w 2 — W3 — W4) T(e — w\ — w 2 + W3 + W4) . 



(17) 



This time the Mellin-Barnes integral does indeed generate poles in e. We choose [32] 



7 



19 



13 



25 



5 5 



(18) 




C2 = "24 



C3 = "24' 



C4 = "48 




in order to separate left poles of T-functions from right ones, and subsequently perform the 
analytic continuation to e = [32]. This generates four kernels, one four-dimensional one, two 
three-dimensional ones, and one two-dimensional one. We arrive at the final result 



Despite the fact that this integral has one more propagator compared to A-j^ it was much simpler 
to evaluate than the former one, and we did not have to introduce any auxiliary parameters 
in integral or series representations, but could proceed until the end by the same techniques as 
described in A^ 2 . The reason for the simplicity of A% is again that it possesses - contrary to 
A-j^ - an outgoing line with an outer vertex where only three lines meet. 

4 Conclusions and Outlook 

In this letter we have evaluated those master integrals for massless three-loop form factors 
which are necessary for the calculation of the fermionic corrections to the quark form factor. 
We obtained analytical results for all coefficients through transcendentality six in the Riemann 
(^-function, as required to obtain the finite part of the form factor at the three-loop level. For 
the integral Aj^ we could obtain a representation which is valid to all orders in e, in terms 
of hypergeometric functions of unit argument. For the other integrals, we derived multiple 
Mellin-Barnes representations from which we extracted all necessary coefficients order by order 
in e. 

The only missing pieces to complete the set of master integrals for massless three-loop form 
factors are therefore the three diagrams in Figure Q] which have nine propagators. It turns 
out that each of them can be expressed in terms of a sixfold Mellin-Barnes representation [42] 
which gives rise to 0(100) single terms upon performing the analytical continuation to e = by 
means of the package MB [32] . Although the number of single integrals is quite large due to the 
extraction of high poles in e, and the evaluation is not completely automated at certain stages, 
the analytical results for these integrals are within reach [42]. 
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A Some technical details on performing MB integrals 



This appendix is devoted to some details about those terms in diagrams A-j^ and A-j^ which 
require the introduction of auxiliary parameters in addition to the MB variables. Consider the 
term 



C2 + i oo cg+ioo 



dw 2 
~27d 



dw 3 
2ni 



3tt 2 csc 2 (7rw 2 )r(w 2 + l)r(l - w 3 )T(-w 3 )T 3 (w 3 y 
(w 2 + l)T(w 2 + w 3 + 2) 



(20) 



Co— 2 OO Co- I OO 



which appears in the computation of Aj <3 at the stage where two MB integrations are left (the 
first MB integration of this integral can be done by Barnes Lemma and corollaries thereof, 
see [28]). The q are as in Eq. (|10p . and we have converted certain combinations of T- functions 
to Csc. 

We perform the integration over w 3 by closing the contour to the right. After summing all 
residues, one obtains 



cn+t oo 



/ 



dwo 



2iri 

C2 — i oo 



^7T 2 CSC 2 (TTW 2 )[lp i0) (w 2 + 2)] 3 57T 4 CSC 2 (7TU;2)'i/ ; ^( U; 2 + 2 



+ 



+ 



2(w 2 + l) 2 4(u; 2 + l) 2 

3tt 2 csc 2 (ttw 2 )[^(w 2 + 2) + 1 e\[1e^ {0) {w 2 + 2) + jj^(w 2 + 2)] 

2(w 2 + l) 2 

7T 2 CSC 2 (7rW 2 )^ (2) (w2 + 2) 7T 2 CSC 2 (7TW 2 )[27| + h^ E TT 2 + 4C 3 ] 



{w 2 + 1) 



4 {w 2 + If 



37r 2 csc 2 (7™ 2 ) 4 F 3 (l, 1,1,1; 2,2, w 2 +3; 1) 



(21) 



(w 2 + l) 2 (w 2 + 2) 

where is Euler's constant. The integration of all but the last term in the above equation can 
be carried out with the standard technique of the nested sums algorithm [36,37]. In the last 
term we write the hyper geometric function as 

i 



dt 



(w 2 +2)(l-t) W2+1 Li 2 (t) 



t 



(22) 
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The integration over W2 can now be done in this term as well. In the end, the integration over 
the auxiliary parameter t can be carried out by means of the package HPL [38,39]. This yields 
+ 3497r 6 /15120 for the expression in Eq. (|2Up , An alternative approach would be to replace 
in Eq. (gJJJ 

l 

n "' 2+l) -rd+T)I dttW2{1 - t)W3 (23) 



r(u> 2 + w 3 + 2) 







with the benefit of having factorized the integrand in W2 and w 3 . The integrations over W2 and 
w 3 can now be carried out by summation of appropriate residues, followed by integration over t. 

As far as Aj^ is concerned, we consider the finite piece of Eq. (|13p where e is set to zero. 
We perform the integration over W4 and close the contour to the right. The result contains two 
hypergeometric functions which we decompose as follows 

4 F 3 (-U;i, -W2,W3 + 1, -Wi -W2 + W3; -Wi - W2, —Wl +W3 + I, -u> 2 + w 3 + 1 ; 1) = 

W Cdt2 Cdtz . r( - 2 +w * + 1)r( - wi +W3+ 1)r( - wi - W2) 



/o Jo Jo T(-wi)T(-W2)T(-w 1 -W2 + w 3 )T 2 (l + w 3 )T(-w 3 ) 

-\l-t 2 ) w * 
(1 - ht 2 t 3 ) 

aF 3 (wi + l,w 2 + l,w 3 + l,w\ + w 2 + w 3 + 2 - wi + w 2 + 2,wi + w 3 + 2,w 2 + w 3 + 2; 1) 

1 j* f\r+ f\u r(w l +W2 + 2)T(w 1 + w 3 + 2)T(w2+w 3 + 2) 

dt\ \ ai2 I dt 3 



Jo Jo 



T( Wl + l)T(w 2 + l)T( Wl +W2 + w 3 + 2)r 2 (l + w 3 )T(-w 3 ) 



t™ 1 (1 - ti) W3 t^ 2 (l - t 2 ) W3 ^1+^2+^3+1(1 _ t3 yw 3 -i 

x (T^P 1 • (25) 

We then perform the integration over w 3 , followed by the other two MB integrations. In the 
end, we carry out the integrations over t±, t 2 and £3. We make simple variable changes in the 
h-t2-t 3 cube where appropriate. 
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